home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48_1 / trap_dn < prev    next >
Internet Message Format  |  1995-03-31  |  3KB
  1. From: David Nurser <itb625nurser@qut.edu.au>
  2. Subject:  v04i032:  trap_dn - Trapezoidal Rule v1.0, Part01/01
  3. Newsgroups: comp.sources.hp48
  4. Followup-To: comp.sys.hp48
  5. Approved: spell@seq.uncwil.edu
  6.  
  7. Checksum: 2883128221 (verify with brik -cv)
  8. Submitted-by: David Nurser <itb625nurser@qut.edu.au>
  9. Posting-number: Volume 4, Issue 32
  10. Archive-name: trap_dn/part01
  11.  
  12. BEGIN_DOC trap.doc
  13. Trapezoidal rule Program  by david nurser (ITB625NURSER@redgum.qut.edu.au)
  14.  
  15. Works exactly the same as the built in integration (via stack)
  16. (almost)
  17.  
  18. Accepts:
  19.  
  20. 1: Lower Limit
  21. 2: Upper Limit
  22. 3: Equation (variable "X")
  23. 4: number of iterations
  24.  
  25. Produces:
  26.  Output...matrix of working (for exams)
  27.           h, the sum, and the answer
  28.  
  29. END_DOC
  30.  
  31. BEGIN_RPL trap.rpl
  32. %%HP: T(3)A(D)F(.);
  33. \<< 'X' PURGE \-> a b
  34. EQ n
  35.   \<< b a - n / 0 \-> h
  36. SUM
  37.     \<<
  38. [[ 0 ]]
  39. { 'n+1' 3 } RDM 0 n
  40.       FOR r r 3 * 1
  41. + r PUT r 3 * 2 + a
  42. r h * + DUP 'X' STO
  43. PUT r 3 * 3 + EQ
  44. \->NUM DUP
  45.         IF r 0 == r
  46. n == OR
  47.         THEN 2 /
  48.         END SUM +
  49. 'SUM' STO PUT 1
  50.       STEP h "h"
  51. \->TAG SWAP SUM "SUM"
  52. \->TAG SWAP SUM h *
  53. "ANS" \->TAG
  54.     \>>
  55.   \>> 'X' PURGE
  56. \>>
  57. END_RPL
  58.  
  59.  
  60. BEGIN_ASC trap.asc
  61. %%HP: T(3)A(D)F(.);
  62. "D9D20E16324563284E20108597632EFE021C432D6E201016D6E201026D6E2020
  63. 5415D6E2010E6E1632D6E201026D6E20101690DA1D6E2010E650FA14B2A21C43
  64. 2D6E201086D6E20303555D4E16328E9209200033920200001000010000000000
  65. 000000000047A208BA20D6E2010E69C2A276BA1B21303F2A2B2130FD0D14B2A2
  66. D6E2010E60A132D6E201027D6E2010273F2A2EEDA19C2A276BA1D6E201027704
  67. D1D6E2010273F2A2EEDA1ED2A276BA1D6E201016D6E201027D6E201086EEDA17
  68. 6BA178BF14563284E20108597632DCC02704D1D6E2010273F2A2EEDA13F2A276
  69. BA1D6E202054154E5A178BF13CE22D6E2010274B2A2279E1D6E201027D6E2010
  70. E6279E1908E1AFE22D9D20ED2A250FA1B21305DF22D6E20303555D476BA14563
  71. 2D6E20303555D497632DCC02704D19C2A208332D6E201086C2A207000086EB52
  72. 2DBBF1D6E20303555D4C2A20B00003555D4EB522DBBF1D6E20303555D4D6E201
  73. 086EEDA1C2A20B000014E435EB522EF532EF5324563284E20108597632EFE029
  74. 3632B213047DF"
  75. END_ASC
  76.  
  77.  
  78. BYTES: #FD74h 397
  79.  
  80. BEGIN_UU trap.uue
  81. begin 644 trap
  82. M2%!(4#0X+46=+>!A(U0V@N0"`5AY-N+O(,$TTN8"`6%M+A`@UN8"`D51;2X0
  83. MX.9A(VTN$"#6Y@(!80FMT>8"`6X%KT$K*L$TTN8"`6AM+C`P5=7D82/H*9`"
  84. M`#,I(````0`0`````````````'0J@*L";2X0X)8L*F>KL1(#\Z*R$@/?T$$K
  85. M*FTN$.`&&B-M+A`@U^8"`7+SHN+>&LFB<K8:;2X0('=`'6TN$"`W+RKNK>$M
  86. M*F>KT>8"`6%M+A`@U^8"`6CNK7&V&H?[064C2"X0@)5G(\T,<D`=;2X0(#<O
  87. M*NZM,2\J9ZO1Y@("15'DI7&X'\,NTN8"`7*THB*7'FTN$"#7Y@(!;G+ID8`>
  88. M^B[2V0+>HE+P&BLQ4/TB;2XP,%75=+8:5#;2Y@(#4U5->3;2S"`'U)$L*H`S
  89. MTN8"`6@L*G```&B^)=*['VTN,#!5U<2B`@L`,%75Y%LBO?O1Y@(#4U5-;2X0
  90. E@.;>&BPJL```04Y3OB7B7R/^-4)E(T@N$("59R/^#I)C(RLQ`"LQ
  91. `
  92. end
  93. END_UU
  94.